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Abstract 


et al. 1999, Bacmann et al. 2002) and in the cold and 


A study of the chemical structure of the envelopes 

nd a sample protostars is introduced — Physical mod- 


arou 

els for the envelopes derived using ID radiative transfer 
modeling of their dust continuum emission are used as in¬ 
put for Monte Carlo modeling of single-dish line observa¬ 
tions in order to establish the chemical inventories for the 
different sources. CO and HCO + abundances are found 
to be correlated with envelope mass: the sources with the 
most massive envelopes show the lowest abundances, sup¬ 
porting the idea that the depletion of these molecules is 
most efficient in the colder, denser environments. Other 
molecules like CS and HCN do not show a similar trend. 
Deuterated species like DCO + and DCN also show signs 
of correlations with envelope mass and temperature: the 
deuterium fractionation of HCO" 1 " is most prominent for 
the coldest and most massive envelopes in the sample and 
at the same time anti-correlated with the fractionation of 
HCN. This puts constraints on the low temperature gas- 
phase deuterium chemistry. 
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1. Introduction 

Understanding the chemistry of protostellar environments 
is of great importance for addressing topics of star for¬ 
mation. It is directly related to processes regulating star 
formation, for example through the ionisation and the 
chemistry also potentially serves as an evolutionary clock 
against which physical models, e.g. for the protostellar 


dense parts of low-mass protostellar envelopes flJprgensen 
et al. 200^) and enhancement of molecules due to libera- 


collapse, can be compared (see e.g. van Dishoeck & Blake 


tion of ice mantles in regions of higher temperatures simi¬ 
lar to “hot cores” or in regions subject to shocks driven by 


protostellar outflows ( 

Bachiller & Perez Gutierrez 19971, 

Ceccarelli et al. 2000, 

schoier et al. 2002). 


(1998) for an overview). Recent evidence suggests that 


In this paper a preliminary overview of the results 
from a large survey of the chemical properties of a sample 
of low-mass protostars is presented. For further details 


see |l0rgensen et al. (2002)| and Jprgensen et al. (2003; in 
prep.). 

2. Sample, observations and modeling 

The sample of objects in this study comprises 18 pre- 
and protostellar objects, including 11 class 0 objects, 5 
class I objects and 2 pre-stellar cores; the full sample is 


presented in Jprgensen et al. (2002). These objects have 


many of the traditional molecular tracers of density and 
temperature may be subject to significant gradients in 
abundances. Examples are the freeze-out of molecules like 


CO onto dust grains both in pre-stellar cores (Caselli 


been observed in a range of molecular transitions using 
the 15 m James Clerk Maxwell Telescope (JCMT), 20 m 
telescope at Onsala Space Observatory and the IRAM 
Pico Valeta 30 m telescope between 2000 and 2002. 

The physical structure of each envelope is derived us¬ 
ing ID radiative transfer calculations. Power-law density 
profiles are assumed for a range of envelope parameters 
(e.g. density slope and optical thickness) and the cor¬ 
responding temperature profiles calculated using the ra¬ 
diative transfer code, DUSTY fllvezic et al. 199^ ) for a 
given central source of heating. For each model, images 
are constructed and compared to the observed SCUBA 
images (450 and 850 fim) and SEDs over the range from 
60 /mi to 1.3 mm. It is found that each protostellar en¬ 
velope can successfully be modeled under these assump¬ 
tions, whereas the pre-stellar cores can not, which is not 
unexpected since these sources per definition are charac¬ 
terized as not having central sources of heating. 

The derived physical models are subsequently used 
as basis for calculating the molecular excitation and ra¬ 
diative transfer for the molecular lines using the Monte 
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Car o code of Hogerheijde fc van der Tak (2000). The 
fractional abundances (assumed to be constant) are con¬ 
strained through comparisons with the integrated line in¬ 
tensities and it is found that the majority of the lines can 
successfully be accounted for with the presented models 
and assumptions. In Fig. |], an overview of the average 
abundances for the class 0 a nd I objects are com pared to 
the molecular cloud L134N ( [Dickens et al. 200C ) and the 
abundances i n the steady state or pre- collapse phase of 
the models of Bergin & Langer (1997). The abundances 
are, where possible, derived from the minor isotopic lines. 
The main isotopic lines of e.g. CS and HCO + are more 
sensitive to departures from spherical symmetry or the 
interaction of molecular outflows with the envelope ma¬ 
terial. 



As can be seen from Fig. 1 and 2 significant differences ex¬ 


ist between the various molecules relative to the two com¬ 
parison sets of abundances and also between the samples 
of class 0 and I objects for some molecules - in particular 
CO and HCO + . CO is thought to be frozen out onto dust 
grains at low temperatures (T < 30 K ; Jprgensen et al 


2002 ) and CO is indeed most depleted for the sources 
with the most massive envelopes (Fig. ||), i.e., the en¬ 
velopes with a region in which the temperatures are low 
enough and densities high enough that the depletion can 
be effective. The primary formation route of HCO + is 
connected to the presence of CO through the reaction 
CO + Hg —► HCO + + H 2 , as is indeed illustrated by the 
data in Fig. ||. 
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Figure 1. Comparisons between average abundances with 
respect to H 2 for different molecules for objects with enve¬ 
lope masses larger or smaller than 0.5 Mq, and with the 
abundances towards the “C” position in the dark cloud 
L134N (Dickens et al. 2000) and the abundances in the 
stage before the collapse in the chemical models of Bergin 
& Langer (1997). All abundances normalized to the abun¬ 
dances of L134N. (From J0rgensen et al. 2003; in prep.) 


3. Molecular depletion 


0.1 1.0 

Mass / Mq 

Figure 2. CO (upper panel) and CS (lower panel) abun¬ 
dances as functions of mass. (From J0rgensen et al. 2002, 
2003; in prep.) 

Other molecules like CS and HCN do not show a sim¬ 
ilar trend with mass, which at first may seem unexpected 
(lower panel of Fig. |j): CS depletion is expected to occur 
at even lower densities than for CO and CS should stay 
bound to the dust grains to even higher temperatures 
than CO due to its larger binding energy to the dust 
grains. Comparison to the models of Bergin & Langer 
(1997)| shows that the CS abundances observed in the 


protostellar environments are indeed lower by more than 
an order of magnitude, which could indicate that the ob¬ 
served lines only probe the region where CS is depleted. 
A detailed comparison with other sulfur-bearing species 
like SO will be needed to fully address this, and similarly 
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FigureS. HCO + abundance as a function of CO abun- 
dance. Symbols as in Fig. [|. (J0rgensen et al. 2003; in 
prep.) 


for the nitrogen-bearing species that also show interest¬ 
ing features - in particular the strikingly low HNC abun¬ 
dances for the protostars compared to other star-forming 
environments. 


4. Deuterium fractionation 


Another interesting aspect of the chemistry in the en¬ 
velopes around protostars is the degree of deuterium frac¬ 
tionation. Gas-phase deuterium fractionation of molecules 
like HCO + is expe cted to be particularly e fficient at low 
temperatures (e.g. Roberts fe Millar 2000 ) due to small 
zero-point energy differences of reactions like HD + Hlj~ —> 
H 2 D + + H 2 . As illustrated in Fig. f|, the highest degree 
of deuterium fractionation is seen to occur in the pre- 
stellar cores and class 0 objects with low temperatures 
(T < 20 K) and high densities. Interestingly, the HCN 
deuterium fractionation seems to be anti-correlated with 
that of HCO + as shown in the lower panel of Fig. [|: for 
HCN deuterium fractionation through other species like 
CHg" (i.e. forming CH 2 D + ) could become increasingly im¬ 
portant in producing DCN at temperatures T ~ 30 K 
(e.g. Turner 2001). Such a scenario could explain the ob¬ 
served anti-correlation between the deuteration of HCO + 
and HCN. 
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Figure). [DCO + ]/[HCO + ] ratio plotted against mass 
(upper panel) and [DCN]/[HCN] plotted against 
[DCO + l/[HCO + ] ratios (lower panel). Symbols as 
in Fig. [|. (Jprgensen et al. 2003; in prep.) 
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